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Abstract 

A new approach to producing a composite material with negative refraction 
index is demonstrated. It is shown that a photonic structure consisting of two 
dielectric materials, with positive and negative dielectric permittivities, can 
support electromagnetic surface waves which exhibit the unusual electromag- 
netic property of left-handedness (or negative refraction index). Depending 
on the dielectric materials, these surface waves localized at the dielectric in- 
terfaces can be either surface plasmons, or phonons. The detailed geometry 
of the structure determines whether this composite left-handed material is 
isotropic, or anisotropic. 

I. INTRODUCTION 

Electromagnetic (EM) properties of materials can be characterized by two macroscopic 
quantities: dielectric permittivity e and magnetic permeability \i. Propagation properties of 
EM waves in the material is determined by e and \i which regulate the relationship between 
the electric field E and the magnetic field B. Those propagation properties may depend on 
the frequency of the waves since both e{u) and fi(u) are, in general, frequency-dependent. 

For the overwhelming majority of materials both e and \i are positive in the propagation 
frequency band. Therefore, in most materials, as well as in vacuum, the relationship between 
E, H, and propagation wave vector k is given by the right hand rule: k ■ [E x H) > 0. The 
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consequence of this is that the group velocity v g = P/U (where P = c[E x H]/4ti is the 
Poynting flux of the wave and U is the wave energy density) of the propagating wave packet 
points in the same direction as its phase velocity v-^ = cuk/\k\ 2 . 

It was first pointed out by Veselago [1] that wave propagation is also possible in materials 
which have simultaneously negative \x and e. Since in such environments the relationship 
between E, H, and k is given by the left-hand rule, such materials are referred to as the Left- 
Handed Materials (LHMs). Their electromagnetic properties are significantly different from 
those of the right-handed materials because the group and phase velocities of electromagnetic 
waves in LHMs oppose each other: k ■ du/dk < 0. For example, an electromagnetic wave 
incident on an interface between the right and left-handed materials stays on the same side 
of the interface normal [1]. In other words, its refracted angle is negative. This property 
gives rise to another name for the LHMs: materials with negative index of refraction [2]. 
Also, the sign of the Doppler effect is reversed in LHMs: an approaching source appears to 
emit higher frequency waves than the receding one. 

LHMs have recently attracted significant attention because of their promise for develop- 
ing the so called "perfect" lenses [3] and low reflectance surfaces [2] . "Perfect" lens enables 
focusing electromagnetic waves to a spot size much smaller than the wavelength A. Conven- 
tional lenses do not permit the image to be significantly sharper than A/2. 

Since materials with \x < do not naturally occur, LHM has to be artificially constructed. 
An LHM in the microwave frequency band was recently constructed [4] as an array consisting 
of metal rods (which provided e < [5]) and split-ring resonators (which provided /i < [6]). 
Photonic structures have also been known to exhibit negative group velocity due to the 
band folding effect [7]. Below I introduce a new concept of making an LHM by creating 
a photonic structure consisting of dielectric regions with e < separated by thin vacuum 
gaps. Left-handedness of such structures is due to the existence of the surface waves at 
the vacuum/dielectric interfaces. Depending on the nature of the negative e material used, 
those surface waves could be either surface plasmons, or surface phonons. 
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Practical implementation of such photonic structures at a micron scale is encouraged by 
the availability of low-loss dielectrics and semiconductors with negative e, including many 
polar crystals such as SiC, LiTaOs, LiF, and ZnSe. The frequency-dependent dielectric 
permittivity of these crystals, given by the approximate formula [8] e(cu) = e^(u 2 —u 2 L ) / {u 2 — 
cuf,), is negative for uj t < u < uj l . Another example of a medium with e < is the free 
electron gas with e = 1 — u 2 /u 2 , where u p is the plasma frequency Its dielectric permittivity 
turns negative for uj < uj p . 

The major innovation introduced in this paper is the use of the dielectric materials with 
small (of order —1) negative dielectric permittivity in order to construct a left-handed com- 
posite material. This is significantly different from earlier work [4] where a periodic arrange- 
ment of metallic components was assembled to achieve left-handedness in the microwave 
frequency range. For microwave, as well as infrared, frequencies the dielectric permittivity 
of metals is essentially e = — oo. 

The organization of the remainder of the paper is schematically shown in Fig. 1. In 
Section II a single dielectric waveguide (SDW) with e c < dielectric cladding (Fig. 1(a)), is 
shown to exhibit left-handedness due to the existence of the surface waves at the vacuum- 
cladding interface. The effective e e fr and /j, c s of the waveguide are shown to be negative (see 
Fig. 2). Two types of surface waves are considered: surface plasmons at the vacuum/plasma 
interface, and surface phonons at the interface between vacuum and silicon carbide (SiC). 
Due to its remarkable property of having the dielectric permittivity e = — 1 at the wavelength 
(10.6/xm) produced by the conventional CO2 lasers, silicon carbide is also proposed as the 
solid state medium for making a near-field "perfect" lens: the tool for enhanced near-field 
imaging in mid-infrared (see Fig. 4). 

Of course, surface waves can only propagate along (and not across) the waveguide walls. 
Section III consider a stack of dielectric waveguides (Fig. 1(b)), which supports waves prop- 
agating in all directions. Only the waves propagating in a limited range of directions, 
however, exhibit the v g ■ < property (where v g and w p h are the group and phase veloc- 
ities). Therefore, the waveguide stack (WS) is not a proper LHM, and the two-dimensional 
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photonic waveguides shown in Fig. l(c,d) are considered in Sections IV,V as the further re- 
finements of the concept. The square lattice photonic waveguide (SLPW) turns out to be a 
highly anisotropic LHM, with the angle a = Z(V p h, v g ) always satisfying n/2 < a < 3n/2 but 
strongly dependent on the propagation direction (see Fig. 6). SLPW is anisotropic even for 
small wavenumbers. This result is explained using the standard perturbation theory. The 
triangular lattice photonic waveguide (TLPW) is found to be a perfectly isotropic LHM for 
small wavenumbers k x>y <C n/d, where d is the lattice periodicity. Recently [9], it has been 
suggested that the group velocity in LHMs is not aligned with the phase velocity, making 
the perfect lens impossible. It is constructively demonstrated that this is not the case even 
for an artificially constructed LHM shown in Fig. 1(d). The main results are summarized 
in Sec. VI. 



II. PROPAGATION OF SURFACE PLASMONS AND PHONONS IN A SINGLE 

DIELECTRIC WAVEGUIDE 

To illustrate how a negative \i can be mimicked using only dielectrics with negative e, 
consider electromagnetic wave propagation in the horizontal (x) direction in a dielectric 
waveguide shown in Fig. 1(a), with a piecewise constant dielectric constant: e = 1 in the 
vacuum channel (for —b<y<b) and e = e c < inside the cladding (for \y\ > b). Consider a 
confined transverse magnetic (TM) wave with non- vanishing components (E x , E y , H z ), and 
assume, by symmetry, that E y = H z = d y E x = at y = 0. Since we are interested in the 
wave propagation along x, introduce integrated over the transverse direction y quantities 

_ rco _ roo 

E x , y = dy E X;V , H z = dy H z . 
Jo Jo 

From Faraday's and Ampere's laws, assuming that E x (x — > oo) = and integrating by 
parts, obtain, correspondingly, 

dE y iuj ~ dH z iu ~ 

— = -n^H,, — = -e eS E y , (1) 
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where the effective dielectric permittivity and magnetic permeability of a dielectric waveg- 
uide are defined as follows: 

1 f°° , ^ icEJy = 0) 
e cS = -~r / dyeE y , /i eff = 1 + ^ '-. (2) 

The definition of the weight-averaged e is intuitive, and /x e fr is defined so as to eliminate the 
longitudinal component of the electric field E x which does not contribute to the power flow 
along the waveguide. 

Equations (1) yield k 2 c 2 /u 2 = fj, e se e s, necessitating that e c g- and /i e fr be of the same 
sign for a propagating wave. Note that /i c ff is a complicated function of the channel width, 
frequency, and the dielectric constant of the cladding, and could generally be either positive 
or negative. Also, the propagating mode is left-handed only if // e fr < 0, e e fr < 0, necessitating 
that the dielectric constant of the cladding be negative. 

Calculating e c g- and /z e ff requires the exact mode structure in the waveguide. The equation 
for H z is 

d (ldH z \ | d (ldH z \ _u^ H 
dx \e dx J dy \e dy ) c 2 

and for the SDW case the harmonic dependence H z oc expi(k x x — out) is assumed. Intro- 



ducing Xv = \Juj 2 /c 2 — k 2 , Xp — \Jk 2 — e c u 2 /c 2 for k x c < cu, the mode structure is given by 
H z = sin XvUl sin Xvb for < y < b and H z = exp [— x P (y — b)] for b < y < oo. It is a surface 
wave because it is localized near the vacuum-cladding interface. Continuity of E x requires 
the continuity of e~ x d y H z across the vacuum-cladding interface at y — b. The dispersion 
relation oo v. s. k x is found by solving the boundary condition equation Xv = ~ e cXv/ tan^fe 
simultaneously with the equations for Xv and Xv Expressions for the effective permittivity 
and permeability are given by 

' 1 1 — Q,o$Xvb\ ( 1 1 — co$Xvb\ 1 



e c ff 



X P Xvsinxvb J \e c Xp XvSmXvb 



, c Xx (I 1 - cos Xvb , 

Ves = 1 y~- h + : T ' 4 

u 2 smxvb \£cXp Xvsmxvb J 

and similar equations are obtained for the sub-luminal mode with k x c > uj. 



Two cladding materials are examined below: a plasmonic material, and a polar dielectric 
silicon carbide (SiC). Plasmonic materials are characterized by the frequency-dependent 
dielectric permittivity e c = 1 — u^/u^ui + ry), where u p is the plasma frequency, and 7 
is the damping constant. At optical frequencies most metals can be considered plasmonic 
materials. Surface plasmons exist at the interface between a plasmonic material and another 
dielectric (or vacuum) with positive dielectric permittivity. Silicon carbide is a low-loss polar 
dielectric which exhibits Restsrahlen: it's dielectric permittivity is negative for frequencies 
ujt < oj < ul, where ujt = 793cm~ 1 and ojl = 969cm" 1 . 

A. Plasmonic Dielectric Cladding 

The dispersion relation and the corresponding e e fj and yU cff are plotted in Fig. 2(a,b) for 
a SDW with plasma-like cladding (e c (u) = 1 — oo^/uo 2 ) and the gap width 2b = c/uo p . The 
propagating surface mode in a SDW exists at the vacuum/plasma interface - therefore it 
is a surface plasmon. This surface plasmon exhibits left-handedness: it's group velocity 
v g = du/dk < is negative, and so are /i e fr < and e cS < 0. The cutoff at u — 0.9lu p is 
caused by the vanishing of the /x e g. This is very different from the cutoff at uj = uj p in a 
homogeneous plasma-like medium which is caused by the vanishing of e e ff. Note that two 
vacuum/plasma interfaces are required to make a surface plasmon left-handed. 

Why is v g < despite t> p h = uj/k > 0? The total Poynting flux P x = cE y H z /A-n along 
the dielectric waveguide is the sum of the fluxes inside the cladding and in the vacuum 
gap. In the gap, E y and H z are in phase, so P x > 0. Inside the cladding E y reverses its 
sign across the vacuum/cladding interface because e c < 0, and the electric displacement 
vector D y = eE y is continuous. H z is continuous across the interface, resulting in P x < in 
the cladding. For a narrow gap, the integrated Poynting flux is negative, and the mode is 
left-handed. 

We emphasize that not any surface wave is left-handed. Achieving left-handedness re- 
quires that (a) the frequency of the mode lies within the stop-band of the cladding (so that 
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e c < 0), (b) there are two interfaces, and (c) the vacuum gap between the interfaces is small. 
The gap between the two claddings with e c < need not be vacuum - it can be filled with 
another dielectric which has a positive dielectric permittivity to enable the existence of the 
surface waves. The cladding with e c < also needs not be plasma-like. In the next section 
the silicon carbide cladding is considered. 



B. Silicon Carbide Dielectric Cladding 

Silicon carbide (SiC) is a polar dielectric [10] with very interesting photonic proper- 
ties because its dielectric permittivity is of order e ~ — 1 for the wavelengths close to 10 
microns. Therefore, a micron-scale waveguide with SiC cladding supporting left-handed 
surface phonons with the frequency corresponding to the vacuum wavelength of 10/im can 
be envisioned. The frequency dependent dielectric permittivity of SiC is given by [10] 

u 2 L - UJ 2 + i-fU 

'c too 2 9 i • ' V / 

UJjp — UJ Z + l^fOJ 

where eoo = 6.7, ujl = 969cm _1 , ujt = 793cm _1 , and 7 = 4.76cm -1 . 

The dispersion relation and the corresponding e e fj and /i c g- are plotted in Fig. 3(a,b) for 
a SDW with the SiC cladding and the gap width 2b = c/uo L . The small damping constant 
7 « 0.005uj l was neglected. The left-handed surface mode is a surface phonon. Qualitatively, 
there is a similarity between the electromagnetic properties of a left-handed plasmon (shown 
in Fig. 2) and those of a left-handed surface phonon (shown in Fig. 3). Both exhibit negative 
group velocity v g = du/dk < and negative e e ff, f^ e s- 

One essential difference is that the group velocity of the surface phonon is much smaller 
than that of the surface plasmon. This is related to the large energy density of the polaritons 
(coupled phonon polarization and photon waves) in SiC. The energy density is proportional 
to U oc d(uje)/duj which is proportional to w| / e 00 /(a'| / — Uj). Both the large value of 
and the narrow width of the reststrahlen band (between ojt and ujl) contribute to the small 
group velocity of the surface phonons. 



7 



Note that there exists a frequency for which both n e g and e c ff are approximately equal 
to (—1). A composite material with /x eff = e c g- = —1 was shown to be ideal for making a 
"perfect" lens [3] which is capable of imaging objects much smaller than the wavelength of 
light. 

Another type of a near-field lens which does not require /i e s = — 1 (although still requires 
e e ff = —1) was also recently suggested [3,11]. It was concluded that a thin slab of material 
with e = — 1 can produce images of sub-wavelength objects with the resolution which by far 
exceeds that of the conventional lenses or even near-field imaging. Unfortunately, finding 
the appropriate material with e = —1 is challenging. Gaseous plasma was suggested as a 
possible candidate for enhanced near-field imaging in the microwave frequency range [11]. 

Unfortunately, the microwave frequency (or the resolution accomplished with mi- 
crowaves) may be too low for many applications. Another promising idea of using thin metal 
layers [3] for enhanced near- field imaging requires extremely thin (less than 100 Angstroms) 
films and very short wavelengths (0.15/zm for Ag films). The reason I've chosen SiC as the 
exemplary ingredient for making a left-handed material is that it may be interesting from the 
technological standpoint: its dielectric permittivity e = — 1 for the frequency corresponding 
to the vacuum wavelength of A = 2-jic/uq = 10.55//m. This wavelength is produced by the 
widely available CO2 lasers. To illustrate the resolution achievable by the combination of 
a C0 2 laser and a thin film of SiC, consider the image of a narrow (sub-wavelength) slit 
produced by a film of SiC of the width d = \ /8rc = 0.42/im. 

Two-dimensional geometry is assumed, and a thin source at x = is assumed to be 
infinitely extended in the z— direction while having the Gaussian shape in the y— direction. 
Transverse magnetic (TM) mode (with field components H Z) E x , and E y ) is assumed, and 
the prescribed magnetic field at the slit H z = exp (—y 2 /a 2 ) with a = A /167r = 0.21/xm is 
assumed. The thin SiC film is located between x = d/2 and x = 3d/2. 

The two-dimensional distribution of the magnetic field H z (x,y) in the rectangular area 
< x < A /27r and — A /2 < y < A /2 is shown in Fig. 4(a). The profiles of H z {y) in three 
planes: at x = (which represents the original shape of the sub- wavelength slit), at x = 2d 
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(the imaging plane), and x — 3d (the SiC slab width behind the image plane), are plotted in 
Fig. 4(b). The enhanced near-field image at x = 2d is practically indistinguishable from the 
shape of the slit at x — 0, and is of much higher fidelity than the un-enhanced near-field at 
x = 3d. Therefore, this numerical example demonstrates the higher fidelity of the enhanced 
near-field imaging in comparison with the standard one. 

To my knowledge, the idea of using a low-loss polar dielectric with the reststrahlen 
band for enhanced near-field infrared imaging with sub-micron resolution is presented here 
for the first time. Possible applications include nano-lithography using a CO2 laser and 
high-resolution biological and chemical sensors. 

III. WAVEGUIDE STACK 

The SDW was considered merely to elucidate the emergence of the left-handed surface 
waves, and to derive the necessary conditions (a-c) for their existence. The objective of 
this Letter is to construct a photonic structure capable of transmitting infinitely extended 
(planar) left-handed waves in all directions. SDW is not infinitely extended in y— direction, 
and it only enables wave propagation in x— direction. Thus, the waveguide stack (WS) with 
periodicity d in y— direction, shown in Fig. 1(b), is considered next. The elementary cell 
of a WS is assumed to have the piecewise constant dielectric permittivity: e(x) = 1 for 
I a; I < b and e(x) = e c for b < \x\ < d/2. Electromagnetic waves in a WS are characterized 

— * 

by the two numbers: phase shift per unit cell — n < <f) y = k y d < ir and the wavenumber k x . 
Electromagnetic waves in a periodically layered medium have been intensely studied since 
the 70's [12], but without the emphasis on surface waves which can exhibit left-handedness. 

Assuming that e c = 1 — u^/u 2 , u p b/c = 0.5, and u p d/c = 4, the dispersion relation in 
the WS uj = u>(ky, k x ) was numerically calculated by integrating Eq. (3) between y = and 
y = d and requiring that H z {y = d)/H z (y = 0) = expi(fi y . It is found that du/dk x < 0, 
duo/dky > 0. Therefore, waves propagate as left-handed along the x-direction, and as right- 
handed along the x— direction. The angle a between the phase and group velocities in 
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the WS is plotted as the dashed line in Fig. 6 v. s. the angle 9 = cos 1 k x /\k\ between the 

— * 

wavenumber k = k x e x + k y e y and the x— axis. For small 9, v g ■ < (a > n/2), and 
the waves are essentially left-handed. However, for large 9, a < ir/2, and the waves are 
right-handed. Therefore, the waveguide stack is not an LHM. 

IV. ELECTROMAGNETIC PROPERTIES OF SQUARE LATTICE PHOTONIC 

WAVEGUIDE 

To ensure that surface waves can propagate in y— direction, an additional set of ver- 
tical channels is introduced, as shown in Fig. 1(c). Rectangular regions of the dielectric 
are separated from each other by a network of orthogonal vacuum channels which form 
a square lattice with the period d = 4.8c/ 'oj p . The channel widths are b = d/8, and the 
dielectric boundaries are rounded with the radius r b = b. Wave propagation along the 
resulting square lattice photonic waveguide (SLPW) is characterized by two phase shifts 
per cell, <p x = k x d and (fi y = k y d. Equation (3) is solved for eigenvalues using a finite 
elements code FEMLAB [13]. The solution satisfies the following boundary conditions: 
H z (d/2,y)/H z (-d/2,y) = e*** and H z (x,d/2)/H z (x,-d/2) = [14]. The local P(x,y) 

— * — * 

and cell-averaged P av Poynting fluxes are computed for each solution. Since v g \\ P av , the 
angle a = l(k, P av ) between P av and k || v p h is used to classify waves into left- and right- 
handed. 

— * 

Two classes of electromagnetic modes were found for any given wavenumber k. The 
higher-frequency mode is right-handed, and the lower-frequency mode is left-handed. The 
physical difference between these modes is illustrated by Fig. 5, where the local Poynting 
fluxes are shown for k x = rf _1 7r/4, k y = 0. The Poynting flux of the left-handed mode 
{uj = 0.85cOp) shown in Fig. 5(a) is localized near the vacuum/cladding interfaces, and points 
in the negative direction. The mode clearly propagates as a surface wave. In contrast, the 
right-handed mode {uj = 0.95uj p ) propagates by tunneling across the dielectric cladding, 
which is why the Poynting flux is positive and not localized, as shown in Fig. 5(b). The two 
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modes have exactly the same frequency cu(k = 0) = cu c = 0.852w p (i. e. are degenerate) for 

— * 

k — 0. As will be demonstrated below, the mode degeneracy results in the strong anisotropy 
of the propagating modes with small (but finite) wavenumbers k: the mode frequency not 

— * — * 

only depends on k = \k\, but also on its direction n — kjk. 

The angle a — l(k, P av ) is plotted in Fig. 6 as the function of the propagation angle 9 for 
the lower-frequency mode of the SLPW. In an isotropic LHM a = it for any 9. In general, 
the medium need not be isotropic, in which case we require that n/2 < a < 3n/2 for all 
propagation directions in order for the medium to be left-handed. 

SLPW is an example of such an anisotropic medium: it is strictly left-handed (because 
7r/2 < a < 3ir/2 for all propagation directions. Group and phase velocities exactly op- 

— # 

pose each other for four different propagation directions: k \\ e x ,e y , (e x ± e y ). Due to the 
anisotropy of the SLPW, radiation "prefers" to propagate along one of the channels. Specifi- 
cally, for k x > k y radiation preferentially flows along x— direction, and vice versa. Symmetry 

— * — * — * 

requires that P || k when k || (e x + e y ), so there is a rapid transition near 9 = n/A from 
the preferential flow along e x for 9 < ir/4 to the preferential flow along e y for 9 > 7r/4. The 
rapidity of the transition is determined by the corner curvature of the structure. Note that 
the anisotropy of the SLPW is not related to the periodicity of the structure. Numerical 
simulations demonstrate that the SLPW is anisotropic even for small \k\ <C d^n far from 
the edges of the Brillouin zone. 

A. Why is SLPW anisotropic for small wavenumbers 

The anisotropy of the SLPW for small k may be somewhat surprising because it implies 

— * 

that the mode frequency is not an analytic function of k. Indeed, if the frequency is a 
quadratic function of k, i. e. u 2 = uj 2 c + gi m kik m (where u c is the cutoff frequency), then the 
symmetry of the square lattice requires that g^i = D5ki, where D is the scalar and Ski is the 
Kronecker delta. This would imply isotropy: from J 2 = oj 2 + Dk 2 follows that the group 
velocity v g = du/dk \\ k is parallel to the phase velocity. However, the frequency is not an 
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analytic function of k. This is due to the fact that at k = the mode is doubly-degenerate 
(see Fig. 5). Below I illustrate using the perturbation theory that the mode degeneracy 
implies anisotropy for the SLPW. 

The solutions of the Eq. (3) corresponding to the phase shift (4> X) <f> y ) per unit cell can 
be expressed as H z = H exp i{k x x + k y y), where if is a periodic function. Substituting H z 
into Eq. (3) yields 

- (vh\ k 2 H . r (2VH ~- A u 2 ~ 

j + — -' k {— + HVe ) = ^ H - (6) 

Equation (6) can be recast as an eigenvalue problem (H^ + + V^)ip = where 
the distances are measured in units of c/uj p , A = uj 2 /u 2 is the eigenvalue, H^ip = — V • 
(e -1 VV0 is the unperturbed Hamiltonian, V^tp = —ik- (2e~ 1 Vtp + tpVe~ 1 ) is the linear in k 
perturbation, and V^ip = k 2 ip/e is the quadratic in A; perturbation of the Hamiltonian. The 
total perturbation V — V± + V 2 . The dielectric permittivity e the function of A. However, for 
simplicity, I will assume in the following that e is frequency independent. For defmitiveness, 
assume that e = —0.4 since, for k = 0, e = 1 — u 2 /u 2 . 

The unperturbed Hamiltonian has a discrete spectrum of eigenfunctions ipu and asso- 
ciated eigenvalues A&. The lowest eigenvalue is A = which corresponds to ipo = const. 
The next doubly-degenerate eigenvalue is = A c = 0.71. The corresponding eigen- 

functions are odd with respect to inversion: ipiji(-x) = —i/)iji(x). They can be conve- 
niently chosen such that ipi(x,—y) = —ipi(x,y), ipj(—x,y) = ipj(x,y), and, correspond- 
ingly, ip n (x,—y) = ipn(x,y), ipn(—x,y) = —ipu(x,y). Note that ipj and ipn are mutually 
orthogonal. The next closest to A c eigenvalue is A 3 = 1.4. This a non-degenerate eigen- 
value, with the corresponding eigenfunction ip 3 satisfying the following symmetry properties: 
fp3(~x, -y) = ip 3 (x,y), ip 3 (x, -y) = -ip 3 (x,y), and i/; 3 (-x,y) = -ip 3 (x,y). 

The ^-dependent perturbation and removes the degeneracy between ipj and ipu, 
resulting in the left-handed and right-handed waves shown in Fig. 5 which have different 
frequencies as indicated earlier. Calculating the mode frequencies in the limit of small k 
can be done in the context of the perturbation theory by solving the secular equation [15]. 
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Defining the overlap integral V^- as 

fd fd 

Vij = dx dy ip*Vipj, 

where ^'s are the unperturbed solutions, the secular equation takes on the form which is 
familiar from quantum mechanics: VkiQ = \ n 5kiC[. The summation over the repeated index is 
assumed, and the sought after n-th perturbed eigenfunction \l/ n is given by \l/ n = J2 m c m>Pm- 
We are primarily interested in the behavior of the perturbed solutions (lower- frequency 

— * 

left-handed mode) and \E' + (upper-frequency right-handed mode) which, in the limit of k — 0, 
become ipi and ipn- 

Strictly, all unperturbed eigenfunctions of H have to be included in the summation. 
For simplicity, I only include the two degenerate eigenfunctions tpi, ipu, and the eigenfunc- 
tions ipo and ip3 which have the closest eigenvalues A = and A = A3. In this truncated 
model, the degenerate doublet of unperturbed eigenfunctions (tpi,tpn) is "surrounded" by 
the two singlets ip and The four eigenfunctions (ip , ipi, ipn, 1P3) are taken to be real and 
normalized to unity. 

Simple symmetry arguments can be used to demonstrate that Vjji = and Vjj = 
Vjiji = Vjj = Ak 2 , where A = (e -1 )/^ is a constant. Thus, the perturbation V does not 
mix if; 1 and ipu, and does not remove the degeneracy of the doublet (ipi,ipu). It is the mixing 
between the doublet and the two singlets via Vi that lifts the degeneracy, and introduces 
the anisotropy of the perturbed solutions and ^/ + . The relevant non-vanishing matrix 
elements are Vj^ = V w , Vjfy = V 20 , 3 = V13, and Vjj^ = V23. It can be shown that 
the effect of the non- vanishing matrix elements Voo = Bk 2 and V33 = Ck 2 (where B and C 

— * 

are constants) on the doublet is of order |/c| 4 . Hence, V 00 and V33 are neglected. Also by 
symmetry, V03 = 0. 

The corresponding secular equation for computing the perturbed frequencies of (\l/ + , , I / _) 
is obtained by solving the following equation [15] for the eigenvalue shift 5\: 
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where A = (<5A) — Ak 2 . In the vicinity of A = 0, the quatric in A Eq. (7) becomes quadratic: 



A 2 + A 



|V in | 2 + |\/ 2 nl 2 . I^l 2 +!^23| 2 \. , |detQ| 2 
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-A c , E 3 = A 3 — A c , and Q is the interaction matrix: 
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Again, symmetry considerations completely determine the functional dependence of the 
matrix elements on k. For example, it can be shown that for the SLPW V10 — iMk y , 
V20 = iMk x , Vi 3 — iNk x , and V 23 = iNk y . Therefore, the bracketed term in Eq. (7) is the 
function of \k\ 2 only. The vanishing of detQ would then guarantee that A is an analytic 
function of k. Not surprisingly, (<5A) = (uj 2 — uj 2 )/uj 2 would be an isotropic function of \k\ 2 . 
For the SLPW, however, | det Q\ 2 = M 2 N 2 {k 2 x — A; 2 ) 2 does not vanish, and 
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(10) 



where the + (— ) signs correspond to ^ + (right-handed) and ^_ (left-handed) modes, re- 
spectively. 

Note from Eq. (10) that A is not an analytic function. It is also manifestly anisotropic 
due to the /c-dependence of Q. The two properties (anisotropy and the non-analyticity) 

— * — * 

go hand in hand: analytic dependence on to on k implies isotropy in the limit of small k. 
In Section V it is found that a periodic lattice with higher symmetry (Triangular Lattice 
Photonic Waveguide) is isotropic. This may be explained by the vanishing of det Q due to 
the higher degree of symmetry of the TLPW. 
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V. ISOTROPY OF TRIANGULAR LATTICE PHOTONIC WAVEGUIDE 



Photonic waveguide can be made much more isotropic if dielectrics are arranged in a 
triangular lattice, as shown in Fig. 1(d), forming a Triangular Lattice Photonic Waveguide 
(TLPW). This arrangement has a higher degree of symmetry than SLPW, and results in 

— * 

a composite medium which is perfectly isotropic for small k well inside the Brillouin zone. 

— * 

The elementary cell of the TLPW and the dispersion relation uo v. s. k for two orthogonal 
directions of k — e XtV k are shown in Fig. 7. The two dispersion curves, which are identical 
for \k\ <C n/d (establishing the isotropy), are drawn to the respective edges of the Brillouin 
zone: < k x d < 2tt/3 and < k y d < 2n/^ [14]. 

These two directions are chosen because they exhibit the maximum possible anisotropy. 

— # 

Indeed, for k = ke x there exists a vacuum/ cladding interface along which the electromagnetic 

— * 

energy can flow in x— direction, Fig. 8(b). No such interface exists for k = ke y , Fig. 8(a). 
Thus, the local Poynting flux patterns are very different for these two directions, as can be 
seen by comparing Figs. 8(a,b). For the parameters of Fig. (8), k — d _1 7r/6, the index of 
refraction is n = —ck/u = —0.2. Index of refraction can be readily tuned in both directions 
by adjusting the parameters of the TLPW: periodicity d and channel width 2b. 

— * 

Despite the differences in the flow patterns, the cell-integrated fluxes P av are identical for 

— * 

all directions of k, and so are the frequencies: u = 0.86a;p. Numerical results unambiguously 
confirm that, for small \k\, phase and group velocities exactly oppose each other. At least 
for this particular LHM, the claim of Valanju et. al. [9] that a = /.(v p h,v g ) ^ n is not 
confirmed. Anisotropy for large \k\ is merely the consequence of the periodicity of the 
photonic structure. 

The angle a = L{k, P av ) between the phase and group velocities was calculated as the 
function of the propagation angle 9 = cos -1 (k x /k). The propagation wavenumber of mag- 
nitude \k\ = <i" 1 7r/4 was rotated by varying < 6 < 2n/3. The angle a for the TLPW 
was found to vary between 0.9947T < a < 1.0067T. This stands in sharp contrast with the 
significant deviation of a from ir for the SLPW as shown in Fig. (6). Therefore, it is con- 
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eluded that, unlike the SLPW, the TLPW is an isotropic (uniaxial) left-handed structure. 
I speculate that the isotropy is related to the fact that det Q given by Eq. (9) vanishes for 
TLPW but not for the SLPW. 

VI. CONCLUSIONS 

In conclusion, several photonic structures supporting left-handed waves were considered. 
The structures consist of dielectrics with the dielectric permittivities e of the opposing signs, 
enabling the propagation of the surface waves. Two types of surface waves, plasmons and 
phonons, were considered. A specific material (SiC) was suggested as the promising candi- 
date for making left-handed photonic structures, as well as for developing a "perfect" lens 
for enhanced near-field imaging. 

The most promising two-dimensional structure with negative index of refraction appears 
to be the Triangular Lattice Photonic Waveguide because it is isotropic (uniaxial) in the wave 
propagation plane. The simpler left-handed structure, the Square Lattice Photonic Waveg- 
uide, was shown to be anisotropic even for infinitesimally small propagation wavenumbers. 
This was shown to be the consequence of the non-analytic dependence of the wave frequency 
on its wavenumber. 

I would like to acknowledge the enlightening conversations with Dr. D. J. Bergman and 
Dr. R. C. McPhedran on the anisotropy of the square lattice photonic waveguide. 
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FIGURES 




FIG. 1. (a) Single dielectric waveguide (SDW) consists of a vacuum gap surrounded by the 
dielectric cladding with e c < 0. (b) One-dimensional waveguide stack (WS). (c) Square lattice 
photonic waveguide (SLPW): dielectric cladding regions (white squares) separated by vacuum 
channels (shaded), (d) Triangular lattice photonic waveguide (TLPW). In (a)-(d): white regions 
represent dielectric cladding, shaded regions represent vacuum channels. 
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FIG. 2. (a) Dispersion relation, and (b) effective dielectric permittivity e e s and magnetic per- 
meability fj, e g. Dielectric permittivity of cladding e c = 1 — ujp/u> 2 , gap width 2b = c/co p . 
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FIG. 3. (a) Dispersion relation, and (b) effective dielectric permittivity e c g and magnetic per- 
meability fi e ff of a dielectric waveguide with SiC cladding. Gap width 2b = c/ujl. 
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FIG. 4. Magnetic field -ff z produced by the Gaussian source of the size a = 0.125c/u;o located 
in the x = plane. A slab of SiC with e = — 1 of width d = 0.25c/uo is located between x = d/2 
and x = 1.5d. The image plane is at x = 2d. (a) Magnetic field H z {x,y) distribution inside the 
simulation domain < x < c/ujo, — Ao/2 < x < Ao/2 (in color online) and (b) lineouts of H z (y) in 
three different planes: x = (source plane, solid line), x = 2d (image plane, dot-dashed line), and 
x = 2.5d (dashed line). 
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FIG. 5. Poynting flux in SLPW for the modes with k x = ci _1 7r/4, k y = 0: (a) lower-frequency 
left-handed mode, P av • k < 0, w/wp = 0.85; (b) higher-frequency right-handed mode, P av ■ k > 0, 
Lu/uip = 0.95. Structure parameters: cell size d = 4.8c/u;p, channel width b = d/4. 
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FIG. 6. Angle between the phase and group velocities a v. s. the propagation angle 
9 = cos^ 1 k x /k, where k = cL~ 1 tt/A. WS parameters: u p d/c = 4, b/d = 1/8. SLPW parame- 
ters: ujpd = 4.8, b/d = 1/8, ri, = b. 
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FIG. 7. Dispersion relation u v. s. \k\ for a triangular arrangement of dielectrics as in Fig. 1(d). 
Equilateral parallelogram with d = ?>c/uj p and (3 = it/3 opening angle - elementary cell of the 
photonic structure. Solid line: k = ke y , < kd < 2^/^/3; dashed line: k = ke x , < kd < 2ir/3. 
Channel widths 2b = 0.6c/u> p , dielectric edges smoothed with radius r& = b. 
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FIG. 8. Power flow in a Triangular Lattice Photonic Waveguide. Structure parameters same 
as in Fig. 7. Wavenumber k = dr l -w/Q, and (a) k = ke y (no vacuum/cladding interfaces parallel 
to k); (b) k = ke x , (interface along k) 
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